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1. INTRODUCTION 

Real power loss reduction is main objective of this work. Various methods conventional [1]-[6] and 
Evolutionary techniques [7]-[19] are applied to solve the problem. Rain Drop Optimization (RDO) Algorithm 
applied to reduce the power loss. Rain drop actions have been imitated to model the algorithm. Natural 
behavior of rain drop is flowing downwards form the peak and it will form small streams during the 
progression from the hill. With reference to the gravitation the raindrop flow as stream then as river from the 
peak of mountains or hill then finally it reaches the sea. Reaching the sea by passing many areas including 
valleys is global optimum. Raindrops are engendered in arbitrary mode in the initial iteration itself, then each 
rain drop will assign a neighborhood by itself then in arbitrary mode neighbor points are produced. Projected 
algorithm begins with arbitrarily engendered solution afterwards exploration has been done sequentially 
around the present point until end of the end criterion and also revision of the present value will be there 
throughout the procedure. Rain Drop Optimization (RDO) Algorithm evaluated in standard IEEE 30, bus test 
system. Voltage deviation and power minimization achieved along with voltage enhancement. 


2. PROBLEM FORMULATION 

Solving the optimal reactive power dispatch (ORPD) problem plays a significant role in the efficient 
operation and planning of the power system. The aim of solving the ORPD is to determine the best operating 
point of system for maximizing the voltage stability, minimizing the system loss and the voltage deviations. 
The best operating point includes the terminal voltages of the generators, taps of transformers and the 
injected reactive powers of the shunt compensators. The solution of the ORPD problem is formulated as an 
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optimization problem applied for assigning set control parameters for a certain objective function, satisfying 
the operating constraints of the system. 


Minimization F(x, y) (1) 


Subject to (2) 


Ey) = 0 (2) 
1% y) =0 (3) 
x = [VG1,.., VGugs QC1,-., Anes Th, --» Tng] (4) 
y = [PGsiacx; Vii se, Vin OG iss QGng; SL rs Sling (5) 


Fitness function (F,) (6) is defined to diminish the power loss (MW), 

F, = Pain = Min [ENT Gm [VŽ + V? — 2 * ViVcosø;j]| (6) 
Fitness function (F,) minimization of Voltage deviation is (7), 

F, = Min [DMP Vix — Vass" + DME Qer — ORE 1] T) 


Fitness function (F3) (8) voltage stability index (L-index) is (9), (10), (11). 


F3 = Min LMax (8) 
Lmax = Maz[L;];j = 1; Nip (9) 
m -INY Fi 1 
And (10) 
Fji =F: =[Yi]* Yal 
Lmax = Max |1 — [Y1] [Y2] x 4 11 
max = Max |1 — yaY] x | an 


Equality constraints 
0 = PG; — PD; — V; Dena Vj [Gycos[®; — Ø;] + Bysin[ø; - ø;]| (12) 
0 = QG; — QD; — V; Djeng V [Gysin[¥, - 9,] + Bycos[D; — 9] (13) 


Inequality constraints 


Pee Paslack S Pysi slack (14) 
og < Qei < Qa at yi S Ng (15) 
VL < VL; < VL™*,i € NL (16) 


Tœ” <T; < 1 i € Ny (17) 
qmin < Qe < Qe, i E€ Nc (18) 


ISLi| < S’ i € Nez (19) 
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VGP < VG; < VGP ,i € Ng (20) 


Multi objective fitness is (21), (22), (23) 


MOF = OF, + x;OF, + yOF; = OF, + [EN x, [VL — VL™]’ + ZNS xg [QG - QGm""]"] + x,0F, (21) 
tes (MOE Vi VL 
VE =) ymin vi, < yrn A 
1 , 1 1 
; GP gaus ggi 
QG” = e a Q L Q A (23) 
QGP", QG; < QG; 


3. RAIN DROP OPTIMIZATION ALGORITHM 

Natural behavior of rain drop is flowing downwards form the peak and it will form small streams 
during the progression from the hill [20]. In the projected algorithm rain drops will be a particle in the 
population, and it has been described as (24), 


DPİ = [Yin Yiv Vers Yit- Yin] i € {1,2,3,.., m} (24) 


6699 


DP? symbolize rain drop, number of variables indicated by “n” and Yik indicate the k" variable in 


(TIEL 


the problem, size of the population is defined by “m”. DP! has been assumed as a vector or point in the “N” 
dimensional axis and through uniform distribution it has been engendered as function with constraints as 
(25). 


Y,x = uniform distribution function (UDF). (Lower,, Upper,) (25) 


With radius “r” around the point of Rain Drop (RD) in the “N” dimensional axis with respect to the 
neighborhood value the change in Rain Drop (RD) value will be modernized. A Rain Drop (RD) “i” which 
possess a neighbor point “j” is indicated by NHP} as (26) and (27), 


||(RD! — NHP;). unit vector,|| < ||radius(r). unit vector,|| (26) 
raduis(r) = raduiSinitia] X f(iteration) (27) 


Naturally the dominant or main (M) neighbor point NHPj, satisfy amongst of all rain drop RD! and 
by objective function (F) the rain drop and neighbor point values are obtained by (28) and (29), 


F(NHPy,) < F(RD') (28) 
F(NHPy,) < F(NHP’) (29) 


No main point for a single rain drop RD! then predominantly in the position of stationary. Then the 
rain drop has been take out of this situation through a procedure of explosion is defined by (30), 


NHPoexploration = neighbor point x explosion base x explosion counter (30) 


The ranking of the rain drops in iteration’s is given by (31) and (32), 


D1i = F(RD')|at titeration — F(RD')|at 1‘titeration (31) 
D2i = F(RD')|at t'"iteration (32) 
order(D1:), order(D2i); are orders at iteration t (33) 


%1, 2 ; are weight cofficients and ranking! is rank of rain drop RD! 


Raindrops are engendered in arbitrary mode in the initial iteration itself, then each rain drop will 
assign a neighborhood by itself then in arbitrary mode neighbor points are produced. Throughout the 
procedure of the optimization the neighbor point will be within in the exploration space limits and if any 
engendered beyond the exploration space limit then it has been modernized by (34) and (35). 
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if(NHP;'), < Lower, then(NHP}), = Lower, (34) 
if(NHP;),. < Upper, then(NHP;),. = Upperr, (35) 


Subsequently for each rain drop and its neighbor points cost function will be computed then the 
comparison of values between rain drop and neighbor point will be done when there is dominant or main 
neighbor point found then the rain drop move towards to that point by altering its position. 

- Begin 

- Parameters values are engendered 

- Fix the each rain drop value as active 

- Rain drop’s initial population produced 

- Iteration =1 

- Drop number =1 

- Explosion counter =0 

- Is rain drop is active? if yes got to next step or else go to step “o” 

- Engender the neighbor points 

- Rain drop and neighbor point cost function value computed 

- Change of rain drop position when a dominant neighbor point has been found then Explosion counter =0 
then go to step “n” 

- Otherwise exploration process will be applied 

- Fix the rain drop as live 

- Swap the rain drop with dominant neighbor point 

- Drop number +1 

- Is drop number <= number of population? If yes go to step “g” 

- Otherwise iteration +1 

- Even after applying the exploration process no neighbor point found means then that particular rain drop is 
marked as motionless 

- After creation of the worth order list then lower order value rain drops are marked as motionless 

- Any live or active rain drop found ; Itertaion <= maximum itertaion 

- If yes go step “fP” 

- Otherwise go to next step “w” 

- For all rain drops the cost function value will be computed 

- Discover the rain drop which possess the minimum cost functional value 

- Then optimal solution is the cost and raindrop position 

- End 


4. SIMULATION RESULTS 

Projected Rain Drop Optimization (RDO) Algorithm evaluated in standard IEEE 30 bus system 
[21]. Active and reactive power consumption is 2.8340 and 1.2620 per unit in 100 MVA base. Table | and 
Table 2 shows the parameters [21]; then Table 3 to Table 4 shows the comparison results. Figures 1-4 gives 
the graphical comparison. 


Table 1. Variables 
System [21] Variables | Minimum (PU) | Maximum (PU) 








VG 0.950 1.10 
Tt 0.90 1.10 
VAR 0.0 5.0 





Table 2. System parameters 








Description [21] IEEE 30 bus 
NB 30 

NG 6 

NT 4 

NQ 9 

NE 41 

Power loss 5.660 

Base casee for VD (Per Uit) 0.582170 
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Table 3. Power loss comparison with different metaheuristicalgorithms 
DE[22] GSA[23] APO-PSOỌ [24] _ RDO 








Power Loss (MW) 4.555 4.5143 4.398 0.940 
Voltage .D (PU) 1.9589 0.87522 1.047 4.226 
L-index (PU) 0.55130 0.141090 0.12670 1.044 





Table 4. Comparison with reference to voltage stability improvement 
DE[22] GSA[23] APO-PSOỌ [24] _ RDO 








Power loss(MW) 6.4755 6.9117 5.698 5.429 
Volatge D (PU) 0.0911 0.0676 0.087 0.083 
L-index (PU) 0.143520 0.13490 0.13770 0.1314 





Table 5. Comparison with reference to voltage deviation minimization 
DE[22] | GSA[23] | APO-PSO [24] RDO 








Power loss (MW) 7.0733 4.9752 4.478 4.249 
Volatge. D (PU) 1.419 0.21579 1.8579 1.8205 
L-index (PU) 0.12460 0.136840 0.12270 0.1173 





Table 6. Comparison of values with reference to multi objective formulation 
APO-PSO [24] _ RDO 

















































































































Power Loss (MW) 4.842 4.738 
Volatge .D (PU) 1.009 1.003 
L-index (PU) 0.1192 0.1184 
5 7 
4 6 
3 E DE [22] 5 E DE [22] 
; m GSA [23] p m GSA [23] 
i m APOPSO [24] 2 m APOPSO [24] 
m RDO 1 m RDO 
0 T T 0 T T 
PLoss VD(PU)  L-index PLoss VD (PU) L-index 
(MW) (PU) (MW) (PU) 
Figure 1. Comparison of power loss Figure 2. Comparison with respect to VSI 
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Figure 3. Comparison of loss with reference with Figure 4. Comparison of loss with reference with 
VDM MOF 


5. CONCLUSION 

In this work optimal reactive power dispatch problem has been successfully solved by Rain Drop 
Optimization (RDO) Algorithm. Projected algorithm begins with arbitrarily engendered solution afterwards 
exploration has been done sequentially around the present point until end of the end criterion and also 
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revision of the present value will be there throughout the procedure. Throughout the procedure of the 
optimization the neighbor point will be within in the exploration space limits and if any engendered beyond 
the exploration space limit then it has been modernized. Rain Drop Optimization (RDO) Algorithm evaluated 
in IEEE 30, bus test system. Voltage stability enhanced, power loss reduced with voltage deviation 
minimization. 
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